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Abstract. We report a numerical evidence that the string tension a can be viewed as an order parameter 
of the phase transition, which separates the smooth phase from the crumpled one, in the fluid surface 
model of Helfrich and Polyakov-Kleinert. The model is defined on spherical surfaces with two fixed vertices 
of distance L. The string tension a is calculated by regarding the surface as a string connecting the two 
points. We find that the phase transition strengthens as L is increased, and that a vanishes in the crumpled 
phase and non- vanishes in the smooth phase. 
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1 Introduction 

A considerable number of studies have been conducted 
on the phase structure of Helfrich and Polyakov-Kleinert 
model of membranes [T||2J|3l|^||oJ | ^||^| [8jf9lll0lll 1II1^II13| - The 
triangulated surface model can be classified into two groups 
|14II15) . One is the model of fixed connectivity surfaces, 
and the other the dynamical connectivity surfaces, which 
are called fluid surfaces. Both kinds of surfaces become 
smooth (crumpled) at infinite (zero) bending rigidity b. 
The model on fixed connectivity surfaces has been consid- 
ered to undergo a finite-6 transition between the smooth 
phase and the crumpled phase. A lot of numerical studies 
including those on fluid surfaces so far support this fact 

seems to be no established understanding of phase transi- 
tions in the fluid surface model. 

Ambjorn et. al. have studied a mass gap and a string 
tension of the fluid model |19|. It was reported in 19 that 
the mass gap and the string tension vanish at the critical 
point of the phase transition, which has been considered 
not to be characterized by a divergence of the specific heat. 
The mass gap was extracted by assuming the spherical 
surface as an oblong one-dimensional string with fixed end 
points separated by a distance L. The string tension was 
also computed by assuming a surface as a sheet of area A 
with fixed boundary, and the same results as those of the 
mass gap were obtained. They used the canonical Monte 
Carlo simulations, which are equivalent with the grand 
canonical ones. 

Recent numerical simulations on the fluid surface model 
suggested that the phase transition is characterized by a 



divergence of the specific heat, although the parameter a, 
the coefficient of the co-ordination dependent term, was 
assumed to have arbitrary values [251126] . Therefore, it is 
interesting to see whether the string tension vanishes or 
not at the critical point of the transition of the model with 
arbitrary a. The notation string tension in this paper cor- 
responds not to the string tension in |19| but to the mass 
gap in |19| : we use string tension in place of the mass gap 
and denote it by a hence force. 

From the simulation studies on the fluid model, we 
obtained a numerical evidence that the string tension a 
vanishes in the crumpled phase and non-vanishes in the 
smooth phase j2Hj- The result presented in [2E] implies 
that a can be considered as an order parameter of the 
phase transition. 

The string tension is considered to be a key to under- 
stand the phase structure of the fluid surfaces. Therefore, 
we show in this paper our simulation data including those 
presented in |28j in order to have an insight into further 
investigations on the phase structure of fluid surfaces. 

We comment on why the result of non- vanishing string 
tension could be a relevant one. It is possible to consider 
that the non-vanishing string tension is connected to two 
interesting problems. The first is the problem of quark 
confinement, which is a problem in high-energy physics. 
The linear potential V(L) ~ L assumed between quark 
and anti-quark separated by distance L gives a finite string 
tension, which is compatible with our result of non- vanishing 
string tension. 

The second is a conversion of external forces into an 
internal energy and vice versa in real physical membranes, 
and is a rather practical problem. We can consider that 



2 



Hiroshi Koibuchi: Grand Canonical simulations of string tension in elastic surface model 



the transition depends on the temperature: the surface be- 
comes crumpled at T>T C and smooth at T <T Cl where T c 
is the transition temperature. Then the surface is picked 
up in two points and extended to sufficiently large at 
T>T C in the beginning; this can be done with zero exter- 
nal force because of the zero string tension. Then, lowering 
the temperature to T <T C , we have a finite string tension 
between the two points on the surface. This is a conversion 
of the internal energy S into an external force. Conversely, 
an external force enlarging the surface in the smooth phase 
can be accumulated as an internal energy. This is also easy 
to understand because the free energy of a string can be 
written as (tension) x (length). If the model in this paper 
represents properties in some real membranes, our result 
implies a possibility of such conversions. 



2 Model 

A sphere in R 3 is discretized with piecewise linear trian- 
gles. Every vertex is connected to its neighboring vertices 
by bonds, which are the edges of triangles. Two vertices 
are fixed as the boundary points separated by a distance 
L. 

The Gaussian energy Si and the bending energy 5*2 
are denned by 

Si =J2( X i~ X jf> S 2 = E( 1 - cos ^)> (!) 

(ij) * 

where Yluj) ^ s the sum over bonds (ij), 9{ in S 2 is the angle 
between two triangles sharing the edge i, and JQ(s R 3 ) 
the position of the vertex i. 

The partition function Z is defined by 

N 

z{b, fttt ;i) = EE / II dx * ex P {~ s ( x i T > Ar )l > ( 2 ) 

N T 1=1 

S(X, T, N)=Si + bS 2 -[iN-a^ log*, 

i 

where denotes the sum over all possible triangula- 
tions T, and N the total number of vertices. It should be 
noted that the chemical potential term —/j,N and the co- 
ordination dependent term —a ^ ■ log qi are included in 
the Hamiltonian. The expression S(X, T, N) shows that 
S explicitly depends on the variables X, T and N. The 
coefficient b is the bending rigidity, and fi the chemical 
potential. Z depends on b, fx, a, and L. The surfaces are 
allowed to self-intersect and hence phantom. 

We consider that the phase structure of the model de- 
pends on the choice of the integration measure J}. qfdXi, 
where qi is the co-ordination number of the vertex i |251 
I26j . The co-ordination dependent term in Eq. J2J) comes 
from this integration measure, because Yii qfdXi can also 
be written as Y\ i dXi exp(a log qi). This a is believed 
to be 2a = 3 |2"31l3()ll31| , and hence it is unclear whether a 
can take arbitrary value. On the other hand qf is consid- 
ered as a volume weight of the vertex i in the integration 



dXi. Thus it is possible to extend a to continuous num- 
bers by assuming that the weight takes a suitable value. 
Therefore, it is interesting to see the dependence of a on 
the phase transitions which can be controlled by the pa- 
rameter a. We note that the continuous a assumed in our 
model does not influence 2a = 3 in the model of [2"§ll3"Ul 

ED 

Note also that the constant term — a log 6 can be 
included in S of Eq. (J2J), because S can be written as 
S = Si + bS 2 — n'N — a J2 i log(gi/6), where fi — fj,'— a log 6. 
As a consequence, the total number N of vertices depends 
on /i and a in the grand canonical simulations using S 
that does not include the constant term. If the simulations 
were done by using S that includes the constant term, the 
results must be equivalent with those without the constant 
term because of the relation between and p! described 
above. 

Let us comment on a relation between the value of a 
and that of q max the maximum co-ordination number, and 
consider why the phase transition is sensitive to the value 
of a. The reason why the phase transition is strength- 
ened at negative a is that the co-ordination dependent 
term — Qf^^og* crumples the surface when a < and 
competes with the bending energy term bS 2 smoothing 
the surface. On the contrary, the term —a log qi tends 
to make q such that q ~ 6 when a > 0. Because of the 
fact that qi is constant on triangulated surfaces due to 
the topological constraint, X^i^Sft becomes maximum 
on the surfaces of uniform co-ordination number q. As a 
consequence, negative a make the surface non-uniform in 
q. Therefore, when a becomes negative large, then g max 
increases, and the surface becomes crumpled. While the 
bending energy S 2 makes the surface smooth, the co-ordination 
dependent term with negative a makes the surface crum- 
pled. Thus two competitive forces co-exist when a is neg- 
ative: one is from the bending energy and the other from 
the co-ordination dependent term. 

We expect 

Z(b, /1, a; L) ~ exp(— <jL) (3) 

in the limit L — > 00 Then, by using the scale invari- 
ance of the partition function, we have |15U19) 

, = M, (4) 
L 

where (Si) and (N) are the mean values of Si and N. 

We note that a surface enclosing two fixed vertices is 
not a one-dimensional string, because the perpendicular 
size of the surface increases with N. However L is chosen 
to be sufficiently larger than the perpendicular size, so 
that Eq. © holds. 

The specific heat, which is the fluctuation of S 2 , is 
defined by C S2 = (b 2 /(N))(d 2 log Z/db 2 ), and is calculated 
by using 

C S2 = -^((S 2 -(S 2 )) 2 ). (5) 
The fluctuation of N denoted by Cn can also be given by 
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(6) 



As we will see later, the phase transition of the model is 
characterized by the divergence of Cs 2 and that of Cjv- 



3 Monte Carlo technique 

X is updated so that X' = X+SX , where 8X takes a value 
randomly in a small sphere. The radius Sr of the small 
sphere is chosen to maintain about 50% acceptance rx for 
the A-update. The radius Sr is defined by using a constant 
number e as an input parameter so that Sr = e (I) , where 
(I) is the mean value of bond length computed at every 
250 MCS (Monte Carlo sweeps). It should be noted that 
Sr is almost fixed because (l) is constant and unchanged 
in the equilibrium configurations. 

T is updated by flipping a bond shared by two trian- 
gles. The bonds are labeled by sequential numbers and 
chosen randomly to be flipped. The rate of acceptance rj- 
for the bond flip is uncontrollable, and the value of r t is 
about 30% < r T < 40%. TV-trials for the updates of X 
and iV-trials for T are done consecutively, and these make 
one MCS. 

N is updated by both adsorption and desorption. In 
the desorption, a vertex is randomly chosen, and then a 
bond that is connected to the vertex is randomly chosen 
so that the two vertices at the ends of the bond unite and 
change to a new vertex. In the adsorption, a triangle is 
randomly chosen in the same way that a bond is chosen 
in the desorption, and a new vertex is added to the center 
of the triangle. As a consequence, the Euler number (=2) 
of the surface remains unchanged in the adsorption and 
the desorption. The acceptance rate rjv is uncontrollable 
as rr is, and the value of r^ is about 55% < r^ < 65% 
in our MC. 

In the adsorption of a vertex, the corresponding change 
of the total energy AS = S (new)— S (old) is calculated. The 
adsorption is then accepted with the probability 
Min[l, exp (-AS) /(N + 1)]. In the desorption, 
AS = S (new)— S (old) is calculated by assuming that one 
vertex is removed. The desorption is then accepted with 
the probability Min [1, TV exp (—AS)]. The adsorption and 
the desorption are tried alternately at every 5-MCS. 

We use surfaces of size N ~ 500, N ~ 1000, and N ~ 
1500. The size N depends on both /j, and a which is fixed 
to three values: a = 5.5, a — 0, and a = —5.5. The reason 
for choosing these three values of a, the phase transition of 
the fluid surfaces is sensitive to a as noted in the previous 
section. The values of fj, are chosen so that TV ~ 500, N~ 
1000, and N ~ 1500 for each a. The diameter L (N) of 
the spheres at the start of MC simulations is fixed so that 
J2t li — 3N/2, where k is the length of the bond i. As a 
consequence, L (N) becomes 



L (N) cx VN 



(7) 



We use three kinds of distance L of the boundary points 
for each Lq(N) such that L — X.5Lq(N), L = 2Lq, and 



L = 3Lo(N). The distance is increased from Lq(N) to L 
in the first 5xl0 6 MCS. 

It should be noted that L = 1.5Lq(N), L = 2Lo, and 
L = 3Lq(N) become oo in the thermodynamic limit TV — > 
oo because of Eq.Q. Therefore, a defined by Eq.© can 
be extracted from these values of L at sufficiently large N. 
Thus, the length L in this paper depends on N and hence 
does not strictly correspond to the one in JUj. In fact, the 
value of L in \nj\ is chosen so that t — L/N changes for 
a given N. However, as we will see, the scaling property 
of physical quantities, such as the dependence of a on N, 
obtained in this paper is compatible with those of a on t 
in H2|. 



4 Results 
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Fig. 1. The normalized histogram h(q) of the co-ordination 
number q at (a) a = — 5.5, /i = —1.87, b = 1.22, (b) a = 
-5.5, /Lt = -1.87, b = 1.34, (c) a = 5.5, fx = 15.9, b = 1.82, and 
(d) a = 5.5, fi — 15.9, b — 1.92. The histograms were obtained 
at the final 2 x 10 7 MCS. 



First, we show in Figs. da)-^d) a normalized his- 
togram h(q) of the co-ordination number q, which is ob- 
tained during the final 2 x 10 7 MCS at (a) a = —5.5, fi = 
-1.87, b = 1.22, (b) a = -5.5, fi = -1.87, b = 1.34, (c) 
a = 5.5, n = 15.9, b = 1.82, and (d) a = 5.5, fi = 15.9, 
b = 1.92. The distance between the two vertices is L = 
3Lo, and the total number of vertices becomes N ~ 1500 
in those cases. The surface becomes crumpled in (a) and 
smooth in (b), and there is no phase transition between 
these phases, as we will see below. We note also that the 
surface becomes crumpled in (c) and smooth in (d), and 
there is a first-order transition between these phases on 
the contrary. We see that the histograms shown in (a) 
and (b) are clearly different from those in (c) and (d). 

In order to show the difference more clearly, we plot 
\ogh(q) against q in Figs. GKa)-Etd). We can see no co- 
ordination number of q > 24 in Figs. [21 a) and[2Jb). To 
the contrary, the curves in Figs. Etc) andEfd) indicate that 
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Fig. 2. Plots of log/i(g) against q at (a) a — —5.5, pi = —1.87, 
b = 1.22, (b) a = -5.5, n = -1.87, b = 1.34, (c) a = 5.5, fi = 
15.9, b = 1.82, and (d) a = 5.5, fx = 15.9, b = 1.92, where h{q) 
are those shown in Figs. |Tf a)-Pfd). 



there exist co-ordination numbers of q > 200 and q > 100 
respectively. The curves of h(q) in Figs. Etc) andEfd) in- 
dicate that configurations of large co-ordination numbers 
play some non-trivial role in the phase transition of fluid 
surfaces. 

The average vertex number (N) are plotted in Figs. 
Eta)-|Bti): (N) in Figs. 02a), 02^b), andEfc) are respectively 
obtained at a = 5.5, L — 1.5Xq; a — 5.5, L = 2Lq; and a — 
5.5, L = 3L . (N) in Figs.^d),^), and^f) are those at 
a = 0, L=1.5£q; a = 0, L = 2Lo; and a = 0, L = 3Lq. (N) in 
Figs.|3Ig),EJh), and|^i) are those at a = — 5.5, L = 1.5L ; 
a = —5.5, L = 2Lq; and a = —5.5, L — 3Lq. The symbols 
A, and O m the figures correspond to those obtained on 
surfaces of size 7V~1000, and 7V~1500 respectively. 

We find from Figs. 02a) EJi) that (N) is weakly depen- 
dent on b and almost independent of L with fixed fi and 
a. The fluctuation Cn of (N) can change against b due to 
this dependence of (N) on b, and will be presented below. 

The specific heat Cs 2 defined by Eq. J5J are plotted in 
Figs. |Ha)-|2£i) an d are respectively obtained at the same 
conditions for (N) shown in Figs. Eta)-Eti) ■ 

Cs 2 at a = 5.5 shown in Figs. Ufa), EJb), and ^c) 
have peaks at intermediate b, however, the growth of peaks 
with increasing N is almost invisible. On the contrary, we 
clearly see the growing of the peaks of Cs 2 at a = 0, L = 
3L in Fig.Eff), and at a = -5.5 in Figs. Efg), Efh), and 
Efi). These indicate that the phase transition strengthens 
not only with decreasing a but also with increasing L at 
least in the region — 5.5<a<0. 

We comment on the total number of MCS and on the 
thermalization MCS. The convergence speed slows down 
when a decreases, because the maximum co-ordination 
number g max increases with decreasing a. 9.6 x 10 s MCS 
were done at a = —5.5, L = 3Lq, b = 1.86, where Cs 2 
has the peak; 7.6 x 10 8 MCS at a = -5.5, L = 2L , b = 
1.86; and 4x 10 s MCS at a = -5.5, L = 1.5£ , b = 1.86. 
Relatively smaller number of MCS was done at b that are 



<N> 
1500 

1000 

500 

<N> 
1500 

1000 

500 



1500 



1000 



500 



(a) |q=5.5, L=1.5L | 



0:H=-l-87 
A:H=-2.25 



q=5.5^ L=2L(7] 




0:H=-1.87 
A:il=-2.25 



(c) q=5.5, L=3L7~| 




0:M=-1.87 
A:n=-2.25 




( g ) a=-5.5, L=1.5L | 


; (h) a=-5.5. L=2L | 


(i) | a=-5.5, L=3L | 


A-A-AAAA 




ooooooo-c 


0:H=15.9 
A:|i=15.5 


0:(i=15.9 
A:n=15.5 


0:ji=15.9 : 
A:n=15.5 









1.7 



1.9 



1.8 



1.9 



1.9 



Fig. 3. The average vertex number (N) against b at (a) a = 5.5, 
L = 1.5L , (b) q = 5.5, L = 2L , (c) a = 5.5, L = 3L , (d) a = 0, 
L = 1.5L , (e) a = 0, L = 2L , (f) a = 0, L = 3L , (g) a = -5.5, 
L = 1.5L , (h) a = -5.5, L = 2L , and (i) q = -5.5, L = 3L . The 
symbols A, and O correspond to those obtained on surfaces 
of size iV~1000, and iV~1500 respectively. 



distant from the transition point, and at a — 0, a — 5.5. 
The thermalization sweeps was about lxlO 7 - 3xl0 7 on 
surfaces of N ~ 1500 at a = —5.5. Relatively smaller MCS 
for the thermalization were done in other cases. 

In order to see the scaling property of the peak values 
we plot Cg^ x against N in log-log scales in Figs. 
0Ja),0Ib), andEfc) respectively obtained at a = 5.5, a = 0, 
and a = -5.5. We find that q? ax at a = 5.5 in Fig. EJa) 
saturate as N increases. On the contrary, C™ ax at a = in 
Fig- OJb) and those at a = —5.5 in Fig. 02c) clearly scale 
according to 



Cmax 
Sr 



C£ ax oc N v . 



(8) 

From the slope of the straight lines in Figs. 02b) and 02c), 
we have 



v = 0.027 ±0.025 [a = 0, L = 1.5L ] 
v = 0.199 ±0.032 [a = 0, L = 3L ] 



(9) 



and 



v = 0.265±0.025 [a = -5.5, L = 1.5L ] , 

v = 0.822 ±0.182 [a = -5.5, L = 3L ]. (10) 

From the value v = 0.199±0.032 at a = 0,L = 3L in Eq. 
© and that ^ = 0.265±0.025 at a= -5.5, L= 1.5L in Eq. 
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Fig. 4. The specific heat Cs 2 against b at (a) a = 5.5, £ = 1.5Lo, 
(b) a=5.5, L = 2L , (c) a=5.5, L = 3L , (d) a = 0, L = 1.5L , 
(e) a = 0, L = 2L , (f) a = 0, L = 3L , (g) a = -5.5, L = 1.5L , 
(h) a = — 5.5, L = 2Lo, and (i) a = —5.5, L = 3Lo- The symbols 
□ , A, and O correspond to those obtained on surfaces of size 
TV ~ 500, TV ~ 1000, and TV ~ 1500 respectively. 



Ijl 0(1 , we understand that the surfaces undergo continuous 
transitions at those conditions. Moreover, ^ = 0.822±0.182 
in Ea. 1(10(1 indicates that the phase transition is of first 
order. 

The peak values C^ ax of the specific heat Cn, which 
is the fluctuation of N defined by Eq. (JBJ, is plotted in 
Figs. Etd), Et e )i andEff). We find also from these figures 
of C]^ ax that the phase transition occurs at a = and 
cv = — 5.5, and that there is no phase transition at cv = 5.5. 
Thus, we confirm that the phase structure described by 
Cn is compatible with that by Cs 2 ■ 

Figures|SJa)-|SIi) are plots of the string tension a against 
b, which were obtained under the conditions that are ex- 
actly same as those in Figs. Ela)-0]1). The string tension 
a is calculated by Eq. |0J. We see in Figs. Efa), [HJb), and 
EJc) that a increases with increasing L and that a on sur- 
faces of small N is relatively larger than that of larger 
surfaces. It is also understood from Figs. EI a ) an d Et D ) 
that a decreases with increasing b on larger surfaces. This 
indicates that a on smooth surfaces are larger than those 
on crumpled surfaces. These properties of a can be seen 
in those obtained at a — in Figs.Etd) andEte), and also 
seen in those obtained at a = —5.5 in Figs.EJg) and EI h)- 
On the contrary, we find from Fig. Eti) that a rapidly 
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Fig. 5. Log-log plots of the peak values C!? ax of the specific 
heat Gs 2 against (TV) obtained at (a) a = 5.5, (b) a = 0, and 
(c) a = —5.5, and log-log plots of the peak values C™ ax of Cn 
against (TV) at (d) a = 5.5, (e) a = 0, and (f) a = —5.5. 



changes at the transition point when TV is increased. We 
already saw in Fig. BJi) that b = 1.86 is the transition 
point of the surface of size TV~1500 at a — — 5.5, L = 3Lq. 
Therefore, we can see in Fig. Eli) that the string tension 
a vanishes in the crumpled phase and non- vanishes in the 
smooth phase. 

In order to see the scaling property of <r, we introduce 
the reduced bending rigidity A such that 



A = 



b c (n,a) 



1. 



(11) 



where b c {/J,,a) denotes the transition point where Cs 2 has 
the peak value shown in Figs. Hta)-|Hi). Then, the tran- 
sition point 6 C ((U, a) is represented by A = 0, the smooth 
phase at b > b c (fi, a) by A > 0, and the crumpled phase 
at b < b c (p , a) by A < 0. The reason why we introduce A 
of Eq. 1(11(1 is because the transition point 6 c (/x, a) moves 
right as TV increases, as confirmed in Fig.^i) for example. 

Figures Ufa)- Efi) show log- log plots of a against (TV) 
obtained at A < 0, A = 0, and A > 0. The straight lines in 
each figure denote the scaling property of a such as 



a (x TV" 



(k > 0). 



(12) 



We confirm from Fig. EJi) that a non- vanishes at A = 0.03 
in the smooth phase, which was expected also from Fig. 
EJi) 32 . Moreover, we find from Figs. Et a ) _ IZt 1 ): an d Eq. 
(113(1 that almost all a satisfy a — ► 0(TV — + oo), which is 
the scaling property at the continuous transition in |19| . 
Recalling that continuous transitions can be seen at a = 
and a = —5.5 in Figs.Etb) and Etc) [orEte) andEtf)], we 
understand that the scaling of a shown in Figs.[7ta)-[7Ii), 
except the non-vanishing <j, are compatible with that in 

Eg. 

The exponent k in Eq. ((12(1 can be obtained by a least 
squares fitting, and some of the results are as follows: 



0.268 ±0.006 [a = 0,L = 3L ,A = 0], 



(13) 



6 



Hiroshi Koibuchi: Grand Canonical simulations of string tension in elastic surface model 



(a) | oc=5.5,L=1.5Lq | 



nn=-3 

A:(i=-2.25 
Q:U=-1.8 



(b) Ot=5.5, L=2l7~| 




(c) ot=5.5, L=3L^| 




1.2 



1.3 



1.2 



1.3 



1.2 



1.3 



(d) | q=Q, L=1.5L I 



□:|i=6 

A:(i=6.7 

Q:|l=7.1 





(f) 0t=0, L=3l7~| 

ooooo 



.1.3 



1.4 1.5 



1.4 1.5 



1.4 1.5 



(g) 



a=-5.5, L=1.5Lq 
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Fig. 6. String tension a against b obtained at (a) a = 5.5, 
L = 1.5L 0) (b) a = 5.5, L = 2L , (c) a = 5.5, L = 3L , (d) a = 0, 
L=1.5L , (e) a = 0, L = 2L , (f) a = 0, L = 3L , (g) a = -5.5, 
L— 1.5io, (h) a = — 5.5, L = 2Lq, and (i) a = —5.5, L = 3Lo- 
The symbols □, A, and O correspond to those obtained on 
surfaces of size /V~500, TV ~ 1000, and TV ~ 1500 respectively. 



k = 0.126 ± 0.019 [a = -5.5, L = 3L , A = -0.03] , 
k = -0.017 ±0.041 [a = -5.5,L = 3L ,A = 0.03] . 



The first k in Eq. 1)13(1 was obtained at a continuous tran- 
sition point, and the second and the third were at the 
discontinuous transition. Although k= —0.017(41) in the 
last of Eq. ((13(1 appears to be ill-defined, we consider that 
it is compatible with the non-vanishing string tension. 

It should be emphasized that the scaling of a in Eq. 
(|T2|) is compatible with a oc {L/N) s with 6 > in [T§] . 
since L oc Lo(N) oc y/N as described in Eq. J7J. Lq(N) is 
the diameter of the initial sphere for the MC simulations 
and chosen to L (N)<x\/N as already noted in Eq. JZJ. 

In fact, we note that S — 2k and S corresponds to 
v/{l — v) in Ref. ^Hj) where v is about v ~ 0.28 in the 
crumpled phase and v ~ 0.42 in the smooth phase close 
to the critical point. These v corresponds to k~0.19 and 
k~0.36 respectively. Thus, these values of n in Ref. ^Hj 
are roughly consistent with the result k — 0.268(6) in Eq. 
1(13(1 . obtained at a continuous transition point. 



(a) | ot=5.5,L=1.5L71 



A :A=-0.03 
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(g) | ot=-5.5, L=1.Il7| 




(h) | ot=-5.5, L=2L~7| 
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□ :X=0.02 





(i) q=-5.5, L=3L7] 
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Fig. 7. Log-log plots of the string tension a against (TV) ob- 
tained at (a) a = 5.5, L = 1.5Lo, (b) a = 5.5, L = 2Lo, (c) 
a = 5.5, L = 3L , (d) a = 0, L = 1.5L , (e) a = 0, L = 2L , (f) 
a = 0, L = 3L , (g) a = -5.5, L = 1.5L , (h) a = -5.5, L = 2L , 
and (i) a =—5.5, L = 3Lo. 



5 Summary and Conclusions 

We have studied the phase structure of the fluid surface 
model of Helfrich and Polyakov-Kleinert by grand canon- 
ical simulations on spherical surfaces with two fixed ver- 
tices of distance L. The model is defined by Hamiltonian S 
containing the Gaussian term Si , the bending energy term 
S2, the co-ordination dependent term S3, and the chemical 
potential term —fj,N: S = Si+frS^— aS^—^N. It is expected 
that the model undergoes a finite-6 transition between the 
smooth phase at b— > 00 and the crumpled phase at b— >0. 
The phase transition was observed at a = and a = —5.5. 
The order of the transition changes from second to first at 
a = — 5.5 with sufficiently large L. The string tension a was 
obtained by regarding the surface as a string connecting 
the two vertices. It is remarkable that a becomes nonzero 
in the smooth phase separated by the discontinuous tran- 
sition from the crumpled phase. Our results indicate that 
a can be viewed as an order parameter of the phase tran- 
sition. It should be noted that our results are compatible 
with those in ^H]) because the obtained a in our study 
vanishes at the critical point of the continuous transition. 

As we have confirmed in this paper, configurations of 
large co-ordination number appear in certain cases and 
play some non-trivial role in the phase transition. Al- 
though we have no clear interpretation of a broad dis- 
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tribution of co-ordination number, it is possible that the 
existence of large co-ordination number is connected with 
some heterogeneous structure of fluid surfaces. 

The results presented in this paper are not conclusive. 
Some problems remain to be studied: Can we find a finite 
string tension in the smooth phase separated by a second- 
order transition from the crumpled one? Can we find that 
the order of the transition remains unchanged on larger 
surfaces? Can we find a clear interpretation of a broad dis- 
tribution of the co-ordination number in biological mem- 
branes? We consider that some points can be resolved by 
the grand canonical MC simulations on sufficiently large 
surfaces. We expect that the non-vanishing string tension 
can also be obtained by the canonical Monte Carlo simu- 
lations on fluid surfaces. Further numerical studies would 
clarify the phase structure of the fluid model of Helfrich 
and Polyakov-Kleinert. 

This work is supported in part by a Grant-in- Aid for 
Scientific Research, No. 15560160. H.K. thanks N.Kusano, 
A.Nidaira, and K.Suzuki for their invaluable help. 
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